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Relic gravitational waves from colliding bubbles and cosmic turbulence 



o 
o 

(N 



> 
oo 

o 
m 
o 
m 
O 

"o ! 

cr 

i ' 



X 



Alberto Nicolis 

Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain 

(Dated: March, 2003) 

A stochastic background of gravitational waves can be generated during a cosmological first order 
phase transition, at least by two distinct mechanisms: collisions of true vacuum bubbles and turbu- 
lence in the cosmic fluid. I compare these two contributions, analyzing their relative importance for 
a generic phase transition. In particular, a first order electroweak phase transition is expected to 
generate a gravitational wave signal peaked at a frequency which today falls just within the band 
of the planned space interferometer LISA. For this transition, I find constraints for the relevant 
parameters in order to produce a signal within the reach of the sensitivity of LISA. The result is 
that the transition must be strongly first order, a > 0.2. In this regime the signal coming from 
turbulence dominates over that from colliding bubbles. 

PACS numbers: 98.80.Cq 



Several phase transitions have probably taken place 
during the evolution of the Universe, and among the pos- 
sible remnants of these cosmic events there is a stochas- 
tic background of gravitational waves (GWs). During a 
first order phase transition the Universe finds itself in a 
metastable phase, the false vacuum, separated from the 
true vacuum by a potential barrier in the free energy of 
an order parameter, usually a scalar field 4>. Quantum 
tunneling throughout the barrier can take place in finite- 
volume regions of space, thus giving rise to the nucleation 
of true vacuum bubbles inside the false vacuum phase. 
As the Universe expands and the cosmic temperature 
T drops down, the energy difference between the false 
and the true vacuum gets larger, and correspondingly the 
phase transition gets more and more probable. Once the 
probability of nucleating one critical bubble per Hubble 
volume per Hubble time becomes 0(1) the transition be- 
gins. This implicitly defines the transition temperature 
T*, by the relation r(T*) ~ H*, where T is the critical 
bubble nucleation rate per unit volume, and is the 
Hubble parameter. The nucleated bubbles begin to ex- 
pand and collide, eventually filling up the whole Universe. 
When two or more bubbles collide spherical symmetry is 
broken, and a small fraction of their kinetic energy is 
released in form of GWs. Moreover, since the bubbles 
expansion is accompanied by macroscopic motions in the 
cosmic plasma, the collision of two (or more) bubbles re- 
sults in some sort of anistropic stirring of the plasma. 
Given the high Reynolds number at such high cosmic 
temperatures, turbulent motions develop, thus providing 
another potentially interesting source of GWs. 

The GW background produced by bubble collisions in 
a first order phase transition has been extensively stud- 
ied in the literature by means of analytical estimates and 
numerical simulations Q> H 0> El Cm the other 
hand, only very recently a systematic computation of the 
amount of GWs emitted by turbulent fluid motions has 
been performed 0, |HJ , showing that turbulence can 
be indeed a powerful source of GWs, even more effective 



than the bubble collision process. The total stochastic 
GW signal one can expect from a generic first order tran- 
sition is the sum of these two contributions: the goal of 
the present paper is to compare them, and to find under 
which conditions one dominates over the other. More- 
over, the electroweak phase transition (EWPT), if it is 
of first order, is expected to give rise to a GW back- 
ground peaked around the milliHertz frequency, and this 
happens to be in the range most relevant for the space in- 
terferometer LISA {n|, which is planned to be launched 
by 2011. I will find what constraints must be satisfied 
by the EWPT in order to give a GW signal within the 
sensitivity of LISA. 

An useful parametrization of a stochastic background 
is given by the dimensionless quantity 
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which is the GW energy density per frequency octave, 
normalized to the critical energy density for closing the 
Universe, p c = 3H§/(8nG), with H = h a x 100 km/(sec- 
Mpc); ho parametrizes the experimental uncertainty in 
Hq. In fact, it is more convenient to characterize the GW 
energy spectrum with ft,Qil gw (/), which is independent on 
ho- The LISA sensitivity to /iQ$7 gw is shown in fig. ^ 

For the purposes of the present paper, we can summa- 
rize the dynamics of a first order phase transition into 
two quantities, traditionally called a and f3. The former, 
a, is defined as the ratio between the false vacuum en- 
ergy (latent heat) density and the plasma thermal energy 
density, both computed at the transition temperature T* . 
It gives a measure of the transition strength: for a <C 1 
the energy jump experienced by the system in performing 
the transition is neglectable with respect to the average 
thermal energy, and thus the transition is very weak; if 
on the contrary a ~ few, the transition is very strongly 
first order. The latter quantity, /3, is the rate of (time) 
variation of the nucleation rate itself, computed at the 
transition time t*, (3 = r/r|*. Therefore /3 _1 fixes the 
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typical time scale of the phase transition: once the tran- 
sition has begun, in a time interval At ~ the whole 
Universe is converted to the true vacuum phase. This is 
because typically (5 is much larger than the Hubble pa- 
rameter at the transition time, (3 3> H* (ljj. In spite of 
the simplicity of their definition, the determination of a 
and (3 for a particular first order phase transition, once a 
specific particle physics model is given, requires involved 
numerical computations (see for instance ref. [lij). 

Two more quantities are important in determining the 
amount of GWs emitted by bubble collisions and turbu- 
lent motions. These are k, the fraction of vacuum energy 
(latent heat) that is transformed into fluid kinetic en- 
ergy, rather than into heat, and Vb, the velocity at which 
the bubble walls expand. For a general phase transition 
their determination involves taking into account frictive 
effects due to departures from thermal equilibrium in the 
vicinity of the wall, which is difficult [ia, llfl ll3 ■ How- 
ever, in the case of a detonation - i.e. if the bubble wall 
propagates faster than the speed of sound (l/v3, in a 
relativistic thermal bath); the opposite case is called de- 
flagration - both Vb and k are functions of a only, in- 
dependently on the particular microphysics that drives 
the transition. Throughout this paper we consider only 
detonations, for the following reasons. First, it has been 
argued that the detonation mode is stable against 
nonspherical hydrodynamical pertubations, while defla- 
gration is not. Second, and more important from our 
perspective, in a deflagration there is no large concentra- 
tion of kinetic energy near the bubble wall, characteristic 
velocities are small, and thus GW production is strongly 
suppressed. We will show that GWs produced at the 
EWPT will be detectable at LISA only if the phase tran- 
sition is very strong, and in such case the bubble expan- 
sion proceeds via detonation Q. For a detonation the 
expansion velocity Vb turns out to be 



v b {a) 



1/V3+ ^Jo? + 2a/3 



1 



(2) 



It is an increasing function of a and varies between the 
speed of sound l/y/3 and the speed of light. Similarly, a 
numerical computation gives [7j 
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which shows that also the efficiency k is increasing with 
increasing a, varying between and 1. 

The energy spectrum of GWs produced by bubble colli- 
sions is peaked at a frequency of about 2/3, independently 
on the bubble expansion velocity Vb • This shows that 
is indeed the relevant time scale of the process. Red- 
shifting that frequency till the present time gives 
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where g* is the number of relativistic degrees of free- 
dom at the transition time. Typical values for (3/H* for 
the EWPT are between 10 2 and a few times 10 3 , with 
T* ~ 100 GeV. This gives a frequency / co n in the range 
(0.5 — 10) mHz, which is precisely the range in which 
LISA achieves its maximum sensitivity. The present time 
fraction of critical density stored in GWs from bubble 
collisions at the peak frequency / co ii is Q 
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At low frequencies the spectrum /iq^coiK/) increases as 
/ 2 ' 8 , while at high frequencies it drops off as J -1 ' 8 0. 

When bubbles collide, the plasma is stirred up at a 
length scale comparable to the size of the colliding bub- 
bles. Larger bubbles are more energetic than smaller 
ones, and indeed it can be shown that the most part of 
the kinetic energy involved in the process is released at 
the largest length scale in the system, namely the radius 
of the largest bubbles at the end of the transition [22J. 
This stirring process creates large eddies in the plasma, 
with a characteristic turnover velocity us given by |l0| 
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Note that, if the phase transition proceeds through det- 
onation bubbles (as we are assuming), also us can be 
expressed uniquely in terms of a, by means of eq. 10. 
Moreover, in such a case it turns out that us is always 
smaller than the bubble wall velocity Vb, eq. 10, and they 
both approach the speed of light in the strong transition 
limit, a ~ few. Once the large eddies have formed, after a 
few revolutions they decay into smaller ones, thus giving 
rise to the usual turbulent energy cascade. The turbu- 
lent motions that develop source a stochastic backgroud 
of GWs peaked at a present (i.e. properly redshifted) 
frequency || 
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1 and thus 



If the transition is strong enough us/vb 

/turb ^ ./coll- 

The relic GW energy density coming from turbulence, 
computed at the peak frequency /turb j is EH 
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At frequencies larger than /t U rb the GW energy spec- 
trum decreases as /~ 7 / 2 , while at lower frequecies it rises 
as f 2 . Strictly speaking, these results holds if magne- 
tohydrodynamical (MHD) effects are absent or negligi- 
ble during the phase transition considered, so that the 
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cosmic fluid develops standard Navier-Stokes turbulent 
motions, with an ordinary Kolmogorov energy spectrum. 
This approach, however, is not always justified: for in- 
stance, large magnetic fields may have been produced 
before or during the EWPT by a number of different pro- 
cesses [13, 0| ■ The inclusion of MHD effects in the above 
GW spectrum is straightforward: one has simply to re- 
place the high frequency behaviour /~ 7 / 2 with /~ 8 / 3 , the 
peak frequency /turb and the peak value of ft.§f^turb being 
unaltered ||. 

Since in a phase transition both bubble collisions and 
turbulent motions may be powerful sources of GWs, it 
is interesting to compare characteristic frequencies and 
energy spectra of the GW backgrounds from the two pro- 
cesses, it order to understand what is the dominant con- 
tribution to a possibly detectable signal. 

First, note that the expected peak frequency of the 
GW background from bubble collisions, eq. 10} , is al- 
ways larger than that of the signal coming from turbu- 
lent motions, eq. (Q), although in the strong transition 
limit ol > 1 they become of the same order. Their ratio 
is indeed 
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if the transition proceeds through detonation bubbles, 
this ratio is a function of a only. 

Then, we compare the relic energy amounts. Given 
their different behaviors with frequency, the dominant 
contribution comes from bubble collisions at high fre- 
quencies and from turbulence at low frequencies. The 
qualitative picture we can expect at intermediate fre- 
quencies is the following, tn moving from lower to larger 
frequencies, we first encounter the turbulence peak at 
/turb; then the spectrum drops down and, when we reach 
/coil, we encounter the bubble-collisions peak. It may 
happen then the second peak is overwhelmed by the high 
frequency tail of the turbulence signal. As we show be- 
low, this possibility depends on the value of a. However, 
even if this is the case, at larger frequencies the contribu- 
tion from bubble-collision is going to become dominant, 
thus giving rise to a change in the slope of the total GW 
spectrum. The opposite case, namely that the turbulence 
peak is overwhelmed by the low frequency tail coming 
from bubble collisions, is in principle possible, but it is 
quantitatively ruled out for any value of the parameters: 
the turbulence peak height is always much larger than 
the height of the low frequency tail from bubble colli- 
sions computed at /turb- Fig. Q schematically shows the 
total GW spectrum for different values of a. 

We can define two interesting quantities. The first is 
the ratio of the GW energy spectra, each computed at 
its own peak frequency, 
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FIG. 1: The total GW spectrum fop^gw vs. frequency (in 
mHz), from a generic first order phase transition, for three 
different values of a, together with the sensitivity curve of 
LISA (bold line) 11]. P/H, is set to 100 and T* to 100 GeV: 
different values would only give a global rescaling of the spec- 
tra, without affecting their shape. 



The second is the same ratio, but with both the spec- 
tra computed at the bubble collision peak frequency; for 
Navier-Stokes turbulence this is 
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while for MHD turbulence it becomes 
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In the case of bubbles expanding as detonation fronts, 
both r\ and depend only on a. In particular, they 
are increasing functions of a, they go to zero in the weak 
transition limit, a — > 0, and get very large in the strong 
transition limit, a £ 1. Therefore, if the phase transi- 
tion is very strongly first-order the bubble-collision peak 
is overwhelmed by the signal coming from turbulence, 
although at higher frequencies it is the bubble-collision 
contribution that dominates the spectrum. This is the 
top case in fig-Q an d happens if T2 > 1, which in terms 
of a reads a > 0.67 in the Navier-Stokes turbulence case 
and a > 0.44 in the MHD one. For smaller values of a 
the spectrum shows both the peaks: the bubble-collision 
one is higher than the turbulence one if r\ < 1, or equiv- 
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alently a < 0.012. This is the bottom case in the figure, 
and corresponds to a weak transition. In the intermedi- 
ate case 0.012 < a < 0.67 (0.44) the turbulence peak is 
higher than the bubble-collision one, see the center plot 
in the figure. 

Finally, we compare the expected GW signals from 
bubble collisions and from turbulence with the LISA sen- 
sitivity curve shown in fig. ^ This easily permits to find 
the values of a and (3 for which the GWs coming from a 
(possibly) first order EWPT are detectable at LISA. To 
do that, we fix T* at 100 GeV [13| and scan the (a,j3) 
plane. For each point considered we superimpose the 
total GW spectrum on the LISA sensitivity curve, and 
check if there is a range of frequencies for which the sig- 
nal is above the instrumental sensitivity and, in such a 
case, which is the background really detectable, either 
the one from turbulence or the one from bubble collisions, 
or both. The results of such a procedure are shown in 
fig. El and described in the caption. The conclusion we 
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FIG. 2: The shaded regions show the values of a and /3/H* 
for which the total GW signal coming from a first order elec- 
troweak phase transition is above the LISA sensitivity curve. 
T» is fixed to 100 GeV. The top plot refers to ordinary Navier 
Stokes turbulence, while the bottom one to MHD turbulence. 
In the region labelled by 1 only GWs from turbulence are 
visible at LISA. In 2 both peaks (from turbulence and from 
bubble collision) are visible. In 3 the bubble collisions peak 
is below the high frequency tail coming from turbulence, but 
nevertheless there exists a frequency range in which the bub- 
ble collision contribution is recognizable by the slope change 
in the signal. In 4 (in the bottom plot) such slope change 
takes place at too high frequecies, where the GW signal is be- 
low the LISA sensitivity, and therefore, the bubble collision 
signal is completely overwhelmed by the turbulence one in the 
frequency window interesting for LISA. 



can draw from the figure is that in order to generate a 
GW signal within the reach of the sensitivty of LISA the 
electroweak phase transition must be strongly first order, 
a > 0.2. In this regime the contribution from turbulence 
to the total GW signal is likely to be more important than 
that from bubble collisions. Such a strong EWPT can 
take place in non-minimal supersymmetric extensions of 
the Standard Model Q- 

I thank Michele Maggiore and Luciano Rezzolla for 
useful discussions. 
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